Magnetocaloric effect in a pyrochlore antiferromagnet Gd2Ti 2 07 



S. S. Sosin, L. A. Prozorova, A. I. Smirnov 
P. L. Kapitza Institute for Physical Problems RAS, 119334 Moscow, Russia 

A. I. Golov, I. B. BerkutcnQ 

Department of Physics and Astronomy, University of Manchester, Manchester M13 9PL, UK 

O. A. Petrenko, G. Balakrishnan 
Department of Physics, University of Warwick, Coventry CV4 7AL, UK 

M. E. Zhitomirsky 

Commissariat a I'Energie Atomique, DSM/DRFMC/SPSMS, 38054 Grenoble, France 

(Dated: October 21, 2004) 

An adiabatic demagnetization process is studied in Gd2Ti207, a geometrically frustrated antifer- 
romagnet on a pyrochlore lattice. In contrast to conventional paramagnetic salts, this compound can 
exhibit a temperature decrease by a factor of ten in the temperature range below the Curie- Weiss 
constant. The most efficient cooling is observed in the field interval between 120 and 60 kOe corre- 
sponding to a crossover between saturated and spin-liquid phases. This phenomenon indicates that 
a considerable part of the magnetic entropy survives in the strongly correlated state. According to 
the theoretical model, this entropy is associated with a macroscopic number of local modes remain- 
ing gapless till the saturation field. Monte Carlo simulations on a classical spin model demonstrate 
good agreement with the experiment. The cooling power of the process is experimentally estimated 
with a view to possible technical applications. The results for Gd2Ti207 are compared to those for 
Gd3GasOi2, a well-known material for low temperature magnetic refrigeration. 



PACS numbers: 75.30.Sg, 75.50.Ee, 75.30.Kz. 

I. INTRODUCTION 

The distinct feature of highly frustrated magnetic ma- 
terials is a peculiar spatial arrangement of the magnetic 
ions. Antiferromagnets on typical geometrically frus- 
trated structures, like, for instance, kagome, garnet, and 
pyrochlore lattices, have an infinite number of classical 
ground states. This macroscopic degeneracy precludes 
any type of conventional magnetic ordering. As a result, 
frustrated magnets remain in a disordered cooperative 
paramagnetic ground state at temperatures well below 
the paramagnetic Curie- Weiss constant 6cw^^ Weak 
residual interactions or quantum and thermal fluctua- 
tions usually induce some kind of ordered or spin-glass 
state at T* <C Oqw- A number of geometrically frus- 
trated magnets have been experimentally studied in the 
past decade^ Magnetic pyrochlore compound Gd 2 Ti 2 07 
is one of the prototype examples. The Gd 3+ ions have 
spin S = 7/2 and zero orbital momentum, which yields 
a good realization of a nearest-neighbor Hciscnbcrg ex- 
change antiferromagnet on a pyrochlore lattice. Re- 
cent specific heat, susceptibility and neutron scattering 
measurements^^ have shown that Gd2Ti2 07 remains 
disordered over a wide temperature interval below its 
6cw — 10 K. The transition to an ordered phase is pre- 
sumably driven by weak dipole-dipole interactions and 
occurs only at T/vi ~ 1 K. 

Infinite degeneracy of the magnetic ground state of a 
frustrated magnet implies the presence of a macroscopic 
number of local zero-energy modes. Such soft modes 



correspond to rotational degrees of freedom of a finite 
number of spins with no change in the total exchange 
energy. In a pyrochlore structure (a matrix of corner- 
sharing tetrahedra on the fee lattice), these groups of 
spins are hexagons formed by the edges of neighbor- 
ing tetrahedra, which lie in the kagome planes [(111) 
and equivalent planes] (see e.g. Ref. 9). In zero ap- 
plied field, if the six spins are arranged antiferromag- 
netically around one hexagon, they effectively decouple 
from the other spins and can rotate by an arbitrary an- 
gle. The low-energy hexagon modes have been observed 
in quasielastic neutron scattering studies on spinel com- 
pound ZnCr2 04, which is a spin-3/2 Heisenberg antifer- 
romagnet on a pyrochlore latticed The thermodynamic 
consequence of these local excitations is that a consider- 
able fraction of the magnetic entropy is not frozen down 
to temperatures much less than 9cw- 

An interesting effect related to a field evolution of the 
zero-energy local modes is an enhanced magnetocaloric 
effect near the saturation field iJ sa t predicted in Ref. OS 
Transformation between a nondegenerate fully polar- 
ized spin state above H sat and an infinitely degenerate 
magnetic state below i? sa t occurs via condensation of a 
macroscopic number of local modes and produces large 
changes of entropy in magnetic field. Quantum fluc- 
tuations do not destroy such an effect In particular, 
at H = i? sat geometrically frustrated magnets on py- 
rochlore and kagome lattices have a finite macroscopic 
entropy, which does not depend on the value of the spin 
of the magnetic ions. The finite entropy of kagome an- 



2 



tiferromagnet has been calculated exactlyjii whereas the 
corresponding quantity for a pyrochlore lattice remains 
unknown. Therefore, an experimental observation of the 
magnetocaloric effect in Gd2Ti2C>7 should be important 
both from a fundamental point of view, as a physical 
probe of the local zero-energy excitations in the frus- 
trated ground state below -ff sa t, and because of the pos- 
sible technological applications. The enhanced magnetic 
cooling power of gadolinium gallium garnet GdsGasO^ 
has been known for a long timejiSai 3 ^ although with- 
out reference to the frustrated nature of a garnet (hyper- 
kagome) lattice. Experiments have been also performed 
on other garnetsJ^ Thus, a comparison between the two 
typical frustrated magnets is important since a stronger 
frustration on a pyrochlore lattice produces the maxi- 
mal cooling rate among all known types of geometrically 
frustrated magnets^ In the present work we study ex- 
perimentally and theoretically the magnetocaloric effect 
in Gd2Ti2C>7. Demagnetization of a Gd2Ti20*7 sample 
under quasi-adiabatic conditions is performed and com- 
pared to a similar study of Gd3Ga50i2- Classical Monte 
Carlo simulations of an ideal adiabatic demagnetization 
process are performed for both materials. Entropy vari- 
ations and the cooling power of Gd2Ti20y in an applied 
magnetic field are calculated from a combination of the 
adiabatic demagnetization results and specific heat mea- 
surements. 



II. EXPERIMENT 

A. Specific heat measurements 

A single- crystal sample of Gd2Ti20y was grown by 
the method described in Ref. [ljj and is approximately 
3.5 x 1.5 x 1 mm 3 in size (32.5 mg by mass). For prelim- 
inary estimation, we present the results of specific heat 
measurements of Gd2Ti207 performed in a Quantum De- 
sign PPMS calorimeter at zero field and at H = 90 kOe 
(see Fig. ^ upper panel) . The curve at H — was taken 
from our previous work in Ref. 18 while the high field 
data were measured in the temperature range from 1.5 
to 20 K. Instead of two sharp peaks on a zero field C(T) 
dependence at Tjvi = 1 K and T/V2 = 0.75 K correspond- 
ing to the ordering phase transitions, the high field curve 
demonstrates a broad maximum around 5-6 K. The lat- 
tice contribution to the specific heat extracted by fitting 
the zero field data at temperatures 30-70 K in the De- 
bye approximation is shown in the figure by a dashed 
line. The temperature dependence of the magnetic en- 
tropy (total entropy with the phonon part subtracted) 
obtained by integrating C/T curves are presented at the 
lower panel of Fig. 1. They give a general overview of the 
magnetocaloric effect in this compound. An adiabatic de- 
magnetization of the system down to zero field starting 
at H = 90 kOe and T = 10 K (around Ocw) results in 
decreasing temperature to 2 K, while the isothermal pro- 
cess at T = 2 K is accompanied by an entropy increase 



by approximately one half of the total magnetic entropy. 
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FIG. 1: Temperature dependence of the specific heat (upper 
panel) and the entropy (lower panel) of the Gd2Ti207 single 
crystal sample: O - zero field, • - H — 90 kOe; solid lines 
are the results of C/T integration, dashed lines represent the 
phonon contribution. 



B. Adiabatic demagnetization 

In the main part of the experiment, the magnetocaloric 
effect in Gd2Ti207 is studied in detail by measuring the 
temperature of the sample in a quasi-adiabatic regime 
as a function of time and magnetic field. For this pur- 
pose, a commercially available thin-film Ru02 resistor 
with the resistance calibrated down to 100 mK in fields 
up to 120 kOe was glued onto the sample. It was also 
used as a heater to regulate the starting temperature of 
the experiment. The sample was suspended on four thin 
constantan wires (20^m in diameter and 5 cm in length) 
soldered to the thermometer to make a 4-wire resistance 
measurement. The experimental cell was put in a vac- 
uum can immersed in a helium bath held at 1.8-4.2 K. 
The heat exchange gas inside the can was absorbed by 
a charcoal cryopump to a pressure of 10~ 7 torr. The 
cryopump was equipped with a heater to desorb some 
exchange gas when necessary to cool the sample during 
the experiment. Magnetic fields up to 120 kOe were gen- 
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FIG. 2: Temperature variations of the adiabatic demagneti- 
zation of Gd2Ti2 07for various starting temperatures; dotted 
lines represent the corrections for the lattice heat capacity 
(see Discussion). Solid lines are obtained by Monte Carlo 
simulations with the exchange constant J = 0.3 K. 

erated by a superconducting magnet. The field H was 
applied perpendicular to the (111) axis. No correction 
for the demagnetization factor was made. The measure- 
ments were done at a sweep rate of 10 kOe/min. 

The temperature versus magnetic field curves for dif- 
ferent starting temperatures are shown in Fig. |5| The 
temperature drop on demagnetization Ti/Tf depends on 
the starting temperature. It has a maximum around 
Ti ~ 10 K where the temperature drops by a factor of 10. 
A characteristic feature for all the Ts(H) curves starting 
below 10 K is that a great deal of cooling occurs in the 
field range from 120 to 60 kOe, which contrasts sharply 
with a continuous adiabatic cooling (T/H = const) of an 
ideal paramagnet. 

Additionally, a parasitic temperature drift T p (due to 
radiation, wiring and the measuring current) was inves- 
tigated at several fixed values of the field for tempera- 
tures along each curve Ts(H). The values of T p have 
a smooth temperature dependence within ±2 mK/sec 
changing sign at about 10 K. The experimental data were 
corrected at each moment of time to (corresponding to 
each value of the magnetic field H) using a set of mea- 
sured T p values by formula 

T CO r(H,to) = T exp (H, to) — / T p dt , (1) 

Jo 

where T exp is the raw experimental value, T cor is the cor- 
rected temperature. One should mention, that the above 
correction procedure is applicable only if the parasitic 
heating does not strongly affect the whole demagneti- 
zation process. In this case, after the correction is ap- 
plied to the raw data, the process appeared to be fully 
reversible by the sweep direction, so that the curves ob- 
tained on increasing the magnetic field do not differ from 



those shown in Fig. [5] A more general heat-balance equa- 
tion 

W p dt = CdT + T {dS/dH) T dH , (2) 

(W p is a parasitic heat leak (T p = W p /C), C is the spe- 
cific heat at a constant field) leads to formula if the 
function T(dS/dH) T /C = -{dT/dH) s does not change 
significantly between T exp and T cor . Such an assumption 
is valid for all the scans shown in Fig. except the one 
starting at Tj = 2 K and the corresponding Ts{H) curves 
are given with the appropriate subtraction. For the low- 
est scan, the correction appears to be of the same order of 
magnitude with the experimental temperature. In addi- 
tion, the derivative (dT/dH)s strongly decreases in the 
vicinity of the ordering transition which leads to an over- 
correction of the experimental data. This overcorrection 
is illustrated in Fig.^by the dashed line. Thus, an ideal 
adiabatic dependence Ts(H) should lie between the raw 
experimental data and the corrected curve. The mini- 
mum temperature reached experimentally on demagneti- 
zation from T t = 2 K at H f = 62 kOe is T min = 0.48 K. 
This cooling limit is associated with the magnetic entropy 
freezing at the transition into an ordered state. When the 
field is further decreased, a weak temperature increase is 
observed. Two temperature plateaus (shown by the ar- 
rows in Fig. 01 are clearly seen in this part of the curve 
corresponding to the phase transitions at Tjvi ~ 1 K and 
T N2 ~ 0.75 K&iLi£. 

It is interesting to compare the above measurements 
to the magnetocaloric effect in another frustrated spin 
system GdaGasO^- The data on Gd3Ga50i2 have been 
obtained on a sintered powder sample (mass 4.05 g) us- 
ing an experimental technique similar to that described 
above. Due to the large heat capacity of the sample, 
no correction for parasitic heat leaks is required in this 
case. The measured adiabatic demagnetization curves for 
Gd3GasOi2 are very similar to the previously published 
results for oriented single crystals^ in the temperature 
range of overlap. Our results are compared to Gd 2 Ti 2 07 
data in Fig. on a field scale normalized to the corre- 
sponding values of H sat - 15 kOe and 70 kOe, respectively. 
The magnetocaloric effect in the Gd3GasOi2 sample is 
qualitatively similar to what is observed in Gd2Ti20y, 
but occurs in a different temperature range due to a 
smaller exchange constant. In a demagnetization pro- 
cess of Gd 3 Ga 5 0i2 starting from H ~ 2_ff sat (28 kOe) 
at T = 0.5 K the temperature decreases by a factor of 
5 reaching its minimum value of 0.1 K in the vicinity 
of the saturation field after which it remains practically 
constant. 



III. DISCUSSION 

Our experimental results demonstrate the decisive 
role of the frustration in producing large adiabatic 
temperature changes in the two materials. In typi- 
cal nonfrustrated antiferromagnets, magnetocaloric ef- 



4 




FIG. 3: Comparison of adiabatic demagnetization in 
Gd2Ti207 (raw data) and Gd3GasOi2 on a field scale nor- 
malized to i/ S at- Monte Carlo simulations are shown by the 
solid lines; the dashed line is an overcorrection to the experi- 
mental data made using formula |T). 



feet in the vicinity of field induced transitions is very 
small. For example, adiabatic AT$ does not exceed 
0.05 K in CsNio.9Feo.1d3i for a similar range of fields 
and temperatures^ For a quantitative description of the 
above data we have performed a series of Monte Carlo 
(MC) simulations on classical antiferromagnetic models 
with only nearest-neighbor exchange. A detailed account 
of the application of the MC technique to model adiabatic 
processes is given in Ref.llOl Quantum spins S = 7/2 are 
replaced by classical n- vectors (|n| = 1), an approxima- 
tion which is valid for high enough temperatures T > JS. 
The only parameter required to compare the MC simula- 
tions to the experimental data is the exchange constant 
J, which is estimated from values of the saturation field 
or the Curie- Weiss temperature. Expressions for these 
quantities for nearest-neighbor antiferromagnets on py- 
rochlore and garnet lattices are 



= 8JS 
= MS 



kB^CW 
kB^CW 



2JS{S- 



1) , 

lJS(S + l), (3) 



5M-B -ffsat 

respectively (see, e.g., Ref. Corresponding val- 

ues for both compounds have been taken from previous 
works^.&iLSi and are summarized in Table HJ For both 
materials, the exchange constants derived from H sat and 
dew agree to within an accuracy of ~ 10 %. The remain- 
ing difference may be attributed to weaker dipolar and 
other interactions, which are responsible for the eventual 
ordering in Gd2Ti207 below Tni and above H = 5 kOe 
in Gd3GasOi2. 

The results of the MC simulations for Gd2Ti207 with 
J = 0.3 K are indicated by the solid lines in Figs. [3 and 
[3J For the scans starting at temperatures above 10 K 
one can not neglect the phonon contribution to the to- 
tal heat capacity of the sample. The influence of lattice 



TABLE I: Experimental values -ff sa t and Qcw and the 
estimations of the exchange constants for Gd2Ti2 07 and 
Gd 3 Ga 5 0i2. 







dew 


./(-Hsat) 


J(8cw) 


Gd 2 Ti 2 7 
Gd 3 Ga 5 0i2 


70 kOe 
15 kOe 


9.5 K 
2.3 K 


0.33 K 
0.1 K 


0.30 K 
0.11 K 



heat on the demagnetization process can be taken into 
account using Eq. J2J, which yields the following rela- 
tionship between the ideal temperature variation AT$ 
and its observed value AT rca i: 



AT S ~ AT rca i(l + Cph/C), 



(4) 



where C p h is the phonon part of the specific heat, C is 
its magnetic part. Using the results shown in Fig. ^ 
one can approximately correct the Ts(H) curves for the 
temperature dependent phonon contribution (see dotted 
lines on Fig. |2J). With this correction, the MC calcu- 
lations demonstrate very good agreement with the ex- 
periment in the whole temperature range at fields above 
-ff sa t. Since a nearest-neighbor Heisenberg pyrochlore an- 
tiferromagnet does not order, the MC simulations cannot 
describe satisfactorily the observed behavior below H sat 
at T < 1 K. Besides, this is also the limit of the valid- 
ity of the classical approximation: JS ~ 1 K. The tem- 
perature increase at H < H sat observed in Gd2Ti207 
(Fig. |3J) may be attributed to the reopening of a gap in 
the excitation spectrum below -ff sa t caused by anisotropic 
interactions. The theoretical dependence obtained for 
Gd3GasOi2 (lower solid line in Fig. [31 with the exchange 
constant J = 0.11 K also fits the high field part of the 
experimental adiabatic curve very well. 

Finally, we estimate the cooling power of the demag- 
netization process using the measured adiabatic Ts(H) 
curves along with the specific heat data C (T) obtained at 
high fields. The amount of heat absorbed by a magnetic 
material during isothermal demagnetization is related to 
the entropy change AQ = TAS\ H f . Consider an adia- 
batic demagnetization curve, which starts at (Hi, Tj) and 
ends at (Hf,Tf). The entropy along the curve Ts(H) re- 
mains constant, so we can relate the entropy variations 
at constant temperature and in constant field: 



AS(T f ) 



Hi 



AS(Hi) 



C(T)/TdT. (5) 



The entropy changes AS of Gd2Ti207 under the 
isothermal demagnetization at various temperatures 
from the initial field Hi = 90 kOe as a function of the fi- 
nal field Hf are presented by closed circles in Fig.0] The 
values of AS are calculated from our experimental data 
using the above equation, while the real entropy of the 
system remains undefined. Two features should be un- 
derlined: (i) about one half of the total magnetic entropy 
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H f (kOe) 

FIG. 4: The molar entropy change of GcbT^Or under isother- 
mal demagnetization from Hi — 90 kOe to Hf at different 
temperatures. Entropy variations of an ideal S — 7/2 para- 
magnet under the same conditions are shown by the dashed 
lines. The inset represents the critical intersection line of the 
two surfaces. 

2R In 8 remains in the system even at temperatures very 
close to the transition into an ordered state at Tjvi = 1 
K; (ii) the largest entropy change AS and consequently, 
a heat absorption AQ occurs in a high field region above 
H sat . This differs significantly from the behavior of an 
ideal paramagnet at low temperatures (shown by grey 
surface on Fig. Q), which releases a significant part of 
its entropy only if demagnetized to Hf <C Hi. This fig- 
ure demonstrates that Gd2Ti2C>7 has a considerable ad- 
vantage below 5 K in the high field cooling power over 
conventional paramagnets. The boundary of the field- 
temperature area where the pyrochlore is advantageous 
over an ideal 5 = 7/2 paramagnet obtained as an in- 
tersection of the two sheets (the bold line) is shown in 
the inset to Fig. For comparison, we performed sim- 
ilar MC simulations of a nearest-neighbor antiferromag- 
net (not shown in the graph). Classical models cannot 
give correct values for the total entropy of a quantum 
spin system. Nevertheless, the variation of entropy is re- 
produced with a remarkable accuracy. For example, the 



calculated values of AS at T = 4.2 K for demagnetization 
from H — 90 kOe to zero field and to Hf = 40 kOe are 
0.52 and 0.40 of the total magnetic entropy respectively. 
The difference from the experiment does not exceed 20%, 
which is due to the large S = 7/2 spin of the Gd 3+ ions. 

The cooling power AQ has a maximum around 4 K 
reaching 30 J/mole Gd. Such an amount of heat corre- 
sponds, for comparison, to the evaporation heat of ap- 
proximately one mole of liquid 3 He at T = 3 K. The 
garnet compound GdsGasO^ is also advantageous as a 
refrigerant material over paramagnetic salts, but at lower 
fields and temperatures. We suggest that a combination 
of the two compounds might become the basis for a inex- 
pensive two-stage adiabatic demagnetization refrigerator 
suitable for effective cooling from T ~ 10 K down to 100 
mK range in a single field sweep. 



IV. CONCLUSIONS 

In conclusion, a large magnetocaloric effect is observed 
in the frustrated pyrochlore antiferromagnet Gd2Ti207 
in agreement with recent theoretical predictions ^fliii 
This observation points at the presence of a macroscopic 
number of local low-energy excitations in the vicinity of 
i/sat. Such modes can be directly probed in quasielas- 
tic neutron scattering measurements analogously to the 
experiment on ZnCr^O^ These excitations may also 
be responsible for the nonfrozen spin dynamics below 
Tjvi observed in Gd2Ti20y by muon spin relaxation 
measurements^ A comparison between our experimen- 
tal data and classical MC simulations shows that this 
numerical technique can semiquantitatively predict the 
magnetocaloric properties of real rare-earth materials 
with large (semiclassical) magnetic moments, which have 
been described so far only in the molecular field approx- 
imation. 
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